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Abstract 

In this work we calculate the corrections to the amputated Green's functions of 4-fermion 
operators, in 1-loop Lattice Perturbation theory. One of the novel aspects of our calculations 
is that they are carried out to second order in the lattice spacing, O (a 2 ). 

We employ the Wilson/clover action for massless fermions (also applicable for the twisted 
mass action in the chiral limit) and a family of Symanzik improved actions for gluons. Our 
calculations have been carried out in a general covariant gauge. Results have been obtained 
for several popular choices of values for the Symanzik coefficients (Plaquette, Tree-level 
Symanzik, Iwasaki, TILW and DBW2 action). 

While our Green's function calculations regard any pointlike 4-fermion operators which 
do not mix with lower dimension ones, we pay particular attention to AF = 2 operators, 
both Parity Conserving and Parity Violating (F stands for flavour: S, C, B). By appropri- 
ately projecting those bare Green's functions we compute the perturbative renormalization 
constants for a complete basis of 4-fermion operators and we study their mixing pattern. 
For some of the actions considered here, even 0(a ) results did not exist in the literature to 
date. The correction terms which we calculate (along with our previous 0(a 2 ) calculation 
of lH3[) are essential ingredients for minimizing the lattice artifacts which are present 
in non-perturbative evaluations of renormalization constants with the RI'-MOM method. 

Our perturbative results, for the matrix elements of AF = 2 operators and for the 
corresponding renormalization matrices, depend on a large number of parameters: coupling 
constant, number of colors, lattice spacing, external momentum, clover parameter, Symanzik 
coefficients, gauge parameter. To make these results most easily accessible to the reader, 
we have included them in the distribution package of this paper, as an ASCII file named: 
4-fermi.m; the file is best perused as Mathematica input. 

The main results of this work have been applied to improve non-perturbative estimates 
of the 5^-parameter in Np = 2 twisted mass lattice QCD [4|]. 

I. INTRODUCTION 

A number of flavour-changing processes are currently under study in Lattice simulations. 
Among the most common examples are the decay K — )■ tttt and K°~K° oscillations. From 
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experimental evidence, we know that these weak processes violate the CP symmetry. In 
theory, the calculation of the amount of CP violation in K°-K° oscillations requires the 
knowledge of the kaon Bk parameter. 

The parameter Bk is obtained from the AS = 2 weak matrix element: 

( ^0|QA5=2|^0 > 

* l(K0\s lti d\0) (0\s llM d\K0) ' {) 

where s and d stand for strange and down quarks, and AS=2 is the effective 4-quark 
interaction renormalized operator, corresponding to the bare operator: 

AS=2 = Hd)(s^dl 7 J = 7,(1 - 75). (2) 

The above operator splits into parity-even and parity-odd parts; in standard notation: 
qAS=2 _ qAS=2^ _0AS=2^__ Since the above weak process is simulated in the framework of 
Lattice QCD, where Parity is a symmetry, the parity-odd part gives no contribution to the 
K°-K° matrix element. Thus, we conclude that Bk can be extracted from the correlator 
(x >0, y <0): 

Ckok(x, y) = ({d75s){x)Ovv+AA( ){dl5s){y)), Oyy+AA = («7a»<*) (s7 M d) + (57^75^(57^75^) , 

(3) 

where Oy V + 2 AA is the bare operator and Oyy+ 2 AA is the corresponding renormalized operator. 

Our results are immediately applicable to other AF = 2 processes of great phenomeno- 
logical interest, such asD — DovB — B mixing. They are also useful in new physics models 
(i.e. beyond the standard model), because there the complete basis of 4-fermion operators 
contributes to neutral meson mixing amplitudes; this is the case for instance of SUSY models 
(see e.g. j|J). For this, one needs to study more general operators of the form 

O xy = (sXd)(sYd) (4) 

where X and Y are general Dirac matrices (see Eq. (jSJ)). 

With Wilson fermions on the lattice, the explicit breaking of chiral symmetry also induces 
a mixing between the Standard Model operator in Eq. (j2j) and the other AS* = 2 operators 
of the basis Eq. (T4j). A strategy which allows to avoid this mixing and, at the same time, 
guarantees automatic (9(a)-improvement of the four fermion operators has been proposed 
in j(| and it makes use of twisted and Osterwalder-Seiler fermions. In this approach, for the 



Bk computation, in place of the operator in Eq. ([3]) a four-quark operator with a different 
flavour content (s, d, s', d'), and with AS 1 = As + As' = 2 is considered, namely 

°vv+aa = (si„d)(s'^d') + (s7 M 7 5 rf)(s / 7^7 5 c/ / ) + (s7/X) {s'l^d) + (&y,sy S (f)(tf'y li 'y B <I) , (5) 

where now the correlator is given by: Ckok'(x,u) = {(d^s)(x)20y V ^ AA (Q)(d'^s')(y)) . 
Making use of Wick's theorem one checks the equality: Ckok'(x,i/) = Ckok( x ,v), which 
means that both correlators contain the same physical information. 

The aforementioned matrix elements are very sensitive to various systematic errors. A 
major issue facing Lattice Gauge Theory, since its early days, has been the reduction of 
effects induced by the finiteness of lattice spacing a, in order to better approach the elusive 
continuum limit. 

In order to obtain reliable non-perturbative estimates of physical quantities it is essen- 
tial to keep under control the 0(a) systematic errors in simulations or, additionally, reduce 
the lattice artifacts in numerical results. Such a reduction, regarding renormalization func- 
tions, can be achieved by subtracting appropriately the 0(a 2 ) perturbative correction terms 
presented in this paper, from corresponding non-perturbative results. 

In this paper we address the perturbative aspects of this problem from a very general 
point of view. In particular, we study the bare 4-point amputated Green's function of the 
most general pointlike 4-fermion operators with four distinct flavours 1 . Although the com- 
putational procedure laid out in the paper is applicable to all orders in the lattice spacing, 
we focus on two different results: 

1. The perturbative 1-loop evaluation of renormalization factors for a variety of 4-fermion 
operators. These factors can be used to renormalize 4-fermion operators computed non- 
perturbatively with any fermion/gluon Wilson-like improved action. This part can be con- 
sidered as an extension of other computations of 4-fermion operator renormalization [7J . 

2. The evaluation of 0(a 2 ) contributions to the aforementioned 1-loop computations. These 
are very useful to improve non-perturbative estimates for the same renormalization fac- 
tors [4j, since they can reduce lattice artifacts, leading to more reliable determinations. 

1 For AS = 1 operators with flavour structure (sXd){qYq) penguin contractions induce a power divergent 
mixing of the four fermion operators with lower dimension operators. This case is not considered in the 
present paper. 
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In Section |TT] we define the general 4-fermion operators and describe the setup of the com- 
putation. The calculations are carried out up to 1-loop in Lattice Perturbation theory and 
up to 0(a 2 ) in lattice spacing. In the same Section we also present simplified expressions for 
the three Feynman diagrams, which constitute the building blocks of the whole calculation. 
In addition, we address certain difficulties which are associated to the 0(a 2 ) computation. 
In Section ILTTI we switch to the evaluation of the renormalization matrices for the 4-fermion 
operators. In particular, we focus on the complete basis of 20 four-fermion operators of 
dimension six which do not need power subtractions (i.e. mixing occurs only with other 
operators of equal dimensions). In the last Section, we summarize the main results of this 
work, and discuss how non-perturbative estimates are being improved by subtracting our 
O (a 2 ) correction terms. 



II. AMPUTATED GREEN'S FUNCTIONS OF 4-FERMION AS = As + As' = 2 
OPERATORS. 

Here we evaluate, up to 0(a 2 ), the 1-loop matrix element of the 4-fermion operators (the 
superscript letter F stands for Fierz.): 

O xy = (sXd){s'Yd') = J2J2 E 4.0*)) (4>)>Wl(*)) (6) 

X c,d hx, k2, k%, &4 

O f xy ee (sXd')(s'Yd) ee Y,Y. E (^)X klk2 d' c k2 (x)) (s> d k3 (x)Y k3k4 d d k4 (x)) (7) 

x c,d k\,k2,k%, ki 

with a generic initial state: d'^(p4) s'" 3 3 (p3)|0), and a generic final state: (0|cZ^ 2 2 (^2) s iHPi)- 
Spin indices are denoted by i, k, and color indices by a, c, d, while X and Y correspond to 
the following set of products of the Dirac matrices: 

X,Y = {l, 1 5 , lfl , lfll 5 ,a llu , 1 5 a tl „}^{S,P,V,A,T,f}; a^ = ^,j„}. (8) 

Our calculations are performed using massless fermions described by the Wilson/clover 
action. By taking m/ = 0, our results are identical also for the twisted mass and the 
Osterwalder-Seiler actions in the chiral limit (in the so called twisted mass basis). For 
gluons we employ a 3-parameter family of Symanzik improved actions, which comprises all 
common gluon actions (Plaquette, tree-level Symanzik, Iwasaki, DBW2, Liischer-Weisz). 
Conventions and notations for the actions, as well as algebraic manipulations involving the 
evaluation of 1-loop Feynman diagrams (up to 0(a 2 )), are described in detail in Ref. 
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di do d3 d4 ds d6 

FIG. 1: 1-loop diagrams contributing to the amputated Green's function of the 4-fermion operator 
O xy . Wavy (solid) lines represent gluons (fermions). 

To establish notation and normalization, let us first write the tree-level expression for the 
amputated Green's functions of the operators Oxy and 0^ Y '- 

A?rL(PliP2,P3,P4, r s ,r d , T a i , T d > )i^Xi4 4 = X hi2 Y hii Ka 2 ^a 3 a 4 , (9) 

(A F )^fa,P2,P3,P4,r a ,r d ,rs,r#)%W^ = -Xi.uY i3 i 2 6 aia J a3a2 , (10) 

where r is the Wilson parameter, one for each flavour. 

We continue with the first quantum corrections. There are twelve 1-loop diagrams that 
enter our 4-fermion calculation, six for each operator Oxy, O^y- The diagrams d\ — d§ 
corresponding to the operator Oxy are illustrated in Fig. [TJ The other six diagrams, d F — d F , 
involved in the Green's function of Oxy are similar to d\ — d$, and may be obtained from 
d\ — dg by interchanging the fermionic fields d and d' along with their momenta, color and 
spin indices, and respective Wilson parameters. 

The only diagrams that need to be calculated from first principles are di, d 2 and d%, 
while the rest can be expressed in terms of the first three. This is a result of the symmetries 
between the diagrams (di, d^), (d 2 , d$) and (d 3 , do). Diagrams d±, d$ and d§ can be expressed 
as diagrams d\, d 2 and d% by exchanging the external quark legs and X, Y, if necessary. In 
particular, the expressions for the amputated Green's functions — can be obtained 
via the following relations: 

A£ y (P1,P2,J*,P4, r„ r d , r sl ,r d ^Z a sU 4 = (A^(-P2. ~Pu -P4, ~Ps, r d , r s , r^.f^Sr)^ 11 ) 
^{PuP2,P3,P4,r s ,r d ,r sl ,r dl )^ l a 2 2 ^ 4 = A^ 2 X (p 3 ,P4,Pi,P2,r s > ,r d > ,r s ,r d )i^£ 2 , (12) 
Ai Y (pi,P2,P3,PA,r s ,r d ,r sl ,r dl )^X!: 4 = ^ 3 X (p 3 ,p 4 , Pl ,p 2 , »V, r rf , r s , r d )%°gg*. (13) 
Once we have constructed A XY — K XY we can use relation: 

(A^fa.M.pa.j^r.^ (14) 
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to derive the expressions for (A F )^ y (i = 1, • • • , 6). From the amputated Green's functions 
for all twelve diagrams we can write down the total 1-loop expressions for the operators 
XY and O f xy : 

Af-^ = £ Aj y , (A F )f-U = E(A F )J y (15) 
j=i i=i 

In our algebraic expressions for the 1-loop amputated Green's functions A^j y , A^ y and 
A^ y we kept the Wilson parameters for each quark field distinct, that is: r s , r d , ry, r> for 
the quark fields s, d, s' and d' respectively. For the required numerical integration of the 
algebraic expressions corresponding to each Feynman diagram, we are forced to choose the 
value for each r parameter. In the numerical results presented in this paper we set: 

r s = r d = TV = r d i = 1. (16) 

Concerning the external momenta pi (shown explicitly in Fig. [T|) we have chosen to 
evaluate the amputated Green's functions at the renormalization point: 

Pl=P2=P3=Pl=P- (17) 

It is easy and not time consuming to repeat the calculations for other choices of Wilson 
parameters and for other renormalization prescriptions. The final 1-loop expressions for 
A^ y , A^ y and A^ y , up to C(a 2 ), are obtained as a function of: the coupling constant g, 
clover parameter csw> number of colors N c , lattice spacing a, external momentum p and 
gauge parameter A. 

As an example we present the results for A^ y and for the special choices: csw = 0, A = 
(Landau Gauge), r s = r d = ry = r# = 1, and tree- level Symanzik improved action: 

Af(p)«^X 4 = tJ^ - * aia ^ 0304 ) x{(A o(a0) )^ + a(A 0(al) )^ +a 2 (A 0(a2) )^| , (18) 



where: 

1 



(Ao(a°))d? — A ilJ2 F l3l4 
+ £(X 7 A Vk* 2 ( y 7'Vk*4 

+ E( x vvk* 2 (*wku 



ln(aV) - 0.05294144(3) 



E(*7 M k , 2 (Y-f)h U [-0.507914049(6)] 



- ln(aV) + 0.018598520(2) 



0.397715726853^^ 
P 



(19) 



The 0(a l ) and 0(a 2 ) contributions of Eq. (fl8|) along with the complete results for all 
diagrams for tree-level Symanzik improved gluons, csw 7^ and A ^ 0, are presented in 



Appendix A; the reader can find similar expressions for other gluon actions in electronic form 
(4-fermi.m). We note in passing that in diagram 3 the dependence on external momentum 
has the same terms as in diagram 2, with identical numerical coefficients; the difference 
between the two diagrams lies in the structure of color and gamma matrices multiplying 
each term. 

The setup presented up to this point applies to both 0(a°) and 0(a 2 ) calculation. For 
the 0(a 2 ) case additional difficulties arise in extracting correctly the full 0(a 2 ) dependence. 
The crucial point of our calculation is the correct extraction of the full 0(a 2 ) dependence 
from loop integrands with strong IR divergences (convergent only beyond 6 dimensions). 
The singularities are isolated using the procedure explained in Ref. In order to reduce 
the number of strong IR divergent integrals, appearing in diagram d\, we have inserted the 
identity below into selected 3-point functions: 

1 = -^{k + ap +k-ap - 2k 2 + 16 V" sm(k a ) 2 sva.(ap a ) 2 ) , (20) 
ap ^ ~' ' 

where q 2 = 4 ^ sin 2 ^) and k (p) is the loop (external) momentum. Repeated use of Eq. 
( 1211]) reduces the 3-point functions to either 2-point functions or more convergent expressions. 
The factor 1/ap 2 in Eq. (|2"U|) can be treated by Taylor expansion. For our calculations it 
was necessary only to O(a ): 

1 - 1 ' ^° P ° +0(a 2 p 2 ). (21) 



ap 2 a 2 p 2 (p 2 ) 2 

Here we present one of the four integrals with strong IR divergences that enter this calcula- 
tion: 

sin (M sin (M - ,f .002457072288 -^¥lU 0.001870841540-^ 



_7r v ; Ar k + ap k — ap F 

\n(a 2 p 2 ) 



a 2 



V( P 2 (o.00055270353(6) 



5127T 2 

ln(a p )\ nnm c^i 0001 n Scr P 



pi (^0.0001282022(1) + -^^r) + 0.000157122310 



P 2 



f , x ln(aV) (pI + pI) 

ViiVA -0.00029731225(4) + ■ y ' - 0.000047949674-^ 



+ 0.000268598599 



768tt 2 ' p 2 

T,gPj 



0(a 4 p)- (22) 



The results for the other three integrals can be found in Ref. [1]. Integrands with simple 
IR divergences (convergent beyond 4 dimensions) can be handled by well-known techniques. 



III. MIXING AND RENORMALIZATION OF Oxy AND O f xy ON THE LATTICE. 



The matrix element (K°\0^= 2 AA \K°) is very sensitive to various systematic errors. The 
main roots of this problem are: a) 0(a) systematic errors due to numerical integration, 
b) with Wilson-like fermions, the operator Oyy^ AA mixes with other 4-fermion AS* = 2 
operators of dimension six. Mixing with operators of lower dimensionality is impossible 
because there is no candidate AS = 2 operator. 

In order to address these problems we have calculated the mixing pattern (renormalization 
matrices) of the Parity Conserving and Parity Violating 4-fermion AS = 2 operators (defined 
below), by using the amputated Green's functions obtained in the previous section. A more 
extensive theoretical background and non-perturbative results, concerning renormalization 
matrices of 4-fermion operators, can be found in Ref. 3 (see also Q, H, Q). Next we 
summarize all important relations from Ref. [s] needed for the present calculation. 

One can construct a complete basis of 20 independent operators which have the symme- 
tries of the generic QCD Wilson lattice action (Parity P, Charge conjugation C, Flavour 
exchange symmetry S=(d <H- d'), Flavour Switching symmetries S'=(s <->• d,s' d') and 
S"=(s <->• d', d ■<->■ s')), with 4 degenerate quarks. This basis can be decomposed into smaller 
independent bases according to the discrete symmetries P, S, CPS', CPS". Following the 
notation of Ref. J&] we have 10 Parity Conserving operators, Q, (P— + 1, S— ± 1) and 10 
Parity Violating operators, Q, (P= — 1, S= ±1): 



of 


=±i 


= \[o vv ±o% 


\+\[o AA ±o F AA ], 




±i = 


\ [O VA ± O f a ] 




Qt 


=±i 


= ±[CW±0£ v 


\-\[o AA ±o F AA ], < 




=±i = 


i [O VA ± o F A ] 


- \ [o AV ± 


Qt 


=±i 


ee \ [Oss ± 0§ s ] 


-\[o PP ±o F PP ], 


Qt 


=±i = 


i [Ops ± O^] 




Qf 


=±i 


ee \ [Oss ± 0§ s ] 


+ ±[Opp±O f p ], 




=±i = 


i [O ps ± O f s ] 




Qt 


=±i 


ee \ [O tt ± Ofr] 


5 


Ql- 




i 

2 







(23) 



Summation over all independent Lorentz indices (if any), of the Dirac matrices, is implied. 
The operators shown above are grouped together according to their mixing pattern. This 
implies that the renormalization matrices Z s=±1 (Z s=±l ), for the Parity Conserving (Vio- 
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lating) operators, have the form: 



rS=±l 



■"11 


Z\2 


"13 


Z14 


zA 

"15 


Z21 


Z22 


Z23 


Z24 


■^25 


Z31 


Z32 


Z33 


Z34 


Z35 


Z41 


Z42 


Z43 


Z44 


Zvo 




Z52 


Z53 


Z54 


Z55 J 



S=±l 



rS=±l 



S=±l 



Z n 

Z22 Z23 

Z32 Z33 

Z 44 £45 

\ 2 M z 55 J 



(24) 



Now the renormalized Parity Conserving (Violating) operators, Q s ±x (Q s ±1 ), are de- 
fined via the equations: 



as=±i _ 7 s=±i n s=±1 n s=±1 — ?s=±\ n s=±1 



(25) 



where l,m = 1,...,5 (a sum over m is implied). The renormalized amputated Green's 
functions L s=±l (C s=±1 ) corresponding to Q s=±1 (Q s=±1 ), are given in terms of their bare 
counterparts L s=±l (C s=±l ) through: 



f S=±l _ y-Z yS=±l rfi=±l fS=±l _ y-Z C2"£>=±1 ft 
u l ~ ^Im ' > *-l — ^* ^Im ' > ~'m 1 



-2r 7 S = ±l T S=±1 rS=±l 



-2~s=±i r S=±i 



(26) 



where Z<q, is the quark field renormalization constant. In order to obtain for a given 
renormalization prescription, one must make use of the inverse fermion propagator, S^ 1 , 
calculated (up to 1-loop and up to 0(a 2 ) for massless Wilson/clover fermions and Symanzik 
improved gluons) in Ref. [lj]. 

The renormalization matrices Z s=±1 (Z s=±1 ), are computed using the appropriate Parity 
Conserving (Violating) Projectors P s=±l (V s=±1 ): 



Pi 



s=±i 



>5=±1 



>S=±1 



pS=±l 



>5=±1 



rw + n 



AA 



QAN C (N C ± 1 



64(iV 2 - 1 



N-aa _j_ n 5 5 — ripp 



32N C (N* - 1) 



n 



vv 



N-aa n 5 5 — ripp 



32iV c (iV 2 - 1) 16(iV 2 - 1) 
II pp 



rigg + ripp 

+ 32jV e (jV e 2 ~l) T 
2JV C ±1 

n S g + ripp 

T 32N C (N? - 1) 



32iV c (iV 2 - 1) 

ripp 

96jV e (iV e 2 -l) ' 
2iV e =Fl 



n 



S=±l 



v s=±i 



Vr 1 



3 

s=±i 



V, 



S=±l 



nvA + n 



64iV c (iV c ±l)' 

riy^ — n^y _^ rigp — Hp5 



64 (^ - 1 



2 -i^ 327V C (7V 2 - 1) 

rivA — n^y rigp — n PS 



32iV c (iV c 2 - 1) 16(iV 2 - 1) 
Tl SP + Tips 



32Jv e (JV|-i) ^ 32NJN? — 1) 

n 5P + rips it 



TT 



= =F 



+ 



TT 



v c y 2Af c Tl 
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where Hxy = (Xizh ® ^443)^0201^0403 • Again, summation is implied over all independent 
Lorentz indices (if any) of the Dirac matrices. The above Projectors are chosen to obey the 
following orthogonality conditions: 



Tr(P, 



s=±i T s=±i 



Tr(7> 



S=±l _ £^=±1 



I ^m{tree)> ~ 11 V I m {tree) / 

where the trace is taken over spin and color indices, and £( tree ) , 



(27) 



" =±1 ^ftretj are the tree-level 
amputated Green's functions of the operators Q s=±1 , Q s=±1 respectively. 



Consistently with the RI' schemes, one may impose the renormalization conditions: 

Tr(P, 5=±1 ■ L* =±1 ) = <U Tr(pf =±1 • t s =^) = 5 lm . (28) 

These conditions should be imposed at a given renormalization scale, /i. Note, however, 
that due to the presence of Lorentz non-invariant quantities, such as Y2 P P% which enter 
the Greens functions at 0(a 2 ) and beyond, the renormalization factors computed in the 
RI'-MOM scheme are also affected through finite cutoff effects by the choice of the direction 
for the external momentum. 

By inserting Eqs. ( 126]) in the above relations, we obtain the renormalization matrices 
Z s=±1 , Z s=±1 in terms of known quantities: 



rS=±l 



7 2 



(D s=±1 ) 



T 



-1 



z 



s=±i 



7 2 



(v s=±i y 



where: 



D 



s=±i 

lm 



Tr(P ; 



S=±l . £ j S=±l\ 



V 



s=±i 

lm 



Tr(P, 



s=±i 



S=±l\ 



Note that D 



s=±i 



and V 



s=±i 



have the same matrix structure as Z 



s=±i 



and Z 



s=±i 



(29) 



(30) 



respec- 



tively. For convenience we express them as: 



D 



s=±i 



16 7T 2 



(df x 


4i 


df 3 


d ± 

"14 


4 


\ 




d ± 

"22 


d ± 

"23 


d ± 

"24 


d ± 

"25 




dfi 


d ± 

"32 


d ± 

"33 


d ± 

"34 






d ± 


"42 


"43 


d ± 


"45 




\4i 


d ± 

"52 


4 


d ± 

"54 


df 5 


J 



+ 0(g 4 



(31) 
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V 



S=±l 



1 + 



16 7T 2 



(*& 





















<*23 


































St 














S M 


si, 


J 



+ 0(g 4 



(32) 



In the parity violating explained in Ref. |8| (Section 5.3), an equality holds 

between two pairs of matrix elements: 



st 2 


= +6 


St 3 


= -5 


sti 


= -5 


si 3 


= +5 



22- 

23; 
32) 
33- 



(33) 
(34) 
(35) 
(36) 



In addition, for the parity conserving projection the matrix elements d£ 3 , d 53 give zero at 
the 1-loop of perturbative theory: 



d 



53 



'53 



0. 
0. 



(37) 
(38) 



The matrix elements of Eqs. fl3T]) - fl32|) have the following simple and generic form: 

,± ... ,±(0,1) , .±(0,2) 2 ,±(0,3) , x ,±(0,4) / ,±(0,5) , ,±(0,6)\ , / 2 2\ 

d l,m - d l,m + C SW % m + C sw d [ m + A Gf, m + (d, m + A G^ m J In (a p J 



2/, ±(2,3) . ,= 

V K,;, + c SW 



±(2,4) 



2 ,±(2,5) . J±(2,6)\ 
+ C SW «Lm + A a Lm ) 



(39) 



+ p m(a p J (d, m + c S w d l m + c sw d, m + A of, m J 



+ 



/i^/ ,±(2,1) 



(C 



Ad 



±(2,2) 
Lm 



) + 0(a 3 ), 



d i,m - d z,m + c sw o l m + c sw d, + A d im + (d im + A d, j m(a p 



+ ^ 



/(C 3) + -swC 4) + 4wC 5) +^C 6) ) 

+ P 2 HaV) (C 7) + ^w C 8) + 4w C' 9) + * 5 

+ ^(C 1) + aC 2) )1 + ^ 3 )- 



(40) 



±(2,10)> 
m , 
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The quantities dff™' and Sf^ appearing above are numerical coefficients depending on the 
number of colors N c and the Symanzik parameters for each gluon action we have considered; 
the index i denotes the power of the lattice spacing a that they multiply. Due to extremely 



lengthy results we provide the quantities d l 



±(i,j) r±(i,j) 

m ' Lm 



(Tables HI- IVIIIj) only for the special 



choices: N c = 3, r s = = r 



■ Td' = 1, and tree- level Symanzik improved action. In all 
Tables the systematic errors in parentheses come from the extrapolation (L — >■ oo) over finite 
lattice sizes. The full set of results is provided in the distribution package of this paper as 
an ASCII file named: 4-fermi.m; the file is best perused as Mathematica input; for notation 
see Appendix B. 

The perturbative renormalization constants (Eqs. ( 1291) ) can be computed directly using 
the 0(a ) coefficients: dfj^ and S^'^ (first line of Eqs. (T4T)1) - ( HIT) ). The renormalization 
factor of the fermion field, Z q , is also required and in the RI'-MOM scheme it reads [3J 



9 2 C F 

16 7T 2 



e« + c sw e (2) + 4 w e^ + Ac« - A ln(«V) 



(41) 



where for the tree-level Symanzik improved action 



,(2) 



r (3) 



-13.0232725(2) 
1.242202721(2) , 
2.01542508(3) , 
4.79200964(9) . 



(42) 
(43) 
(44) 
(45) 



As an example, we provide the exact expression for Z VJ ^\ V) up to 1-loop approximation: 



rS=+l 
'VA+AV 



16 7T 2 



St i 0,1) + csw Sir' + 4 W Str' + A ^| u ' 4; + 2 ln(aV 



r+(0,2) 



(0,3) 



f+(0,4) 



9 2 C, 

16 7T 2 



6«+ Cs we {2) +c| w e( 3 ) + Ae( 4 )-Aln( 



2 2^ 

a p : 



(46) 



This costant is the one relevant for the renormalization of Bk in the twisted 
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mass/Osterwalder-Seiler approach [6], implemented in Ref. j^] 2 . For the tree- level Symanzik 
paction adopted in thecal of Ref. fl, with csw = and in the Landau 
gauge, Eq. (l4"oT) reads: 



VA+AV 



1 - 



16vr 2 



2 ln(aV) +42.3359 



(47) 



While the RI'-MOM scheme allows for a non-perturbative renormalization procedure of 
the lattice operators, the Wilson coefficients entering the effective weak Hamiltonian for 
neutral meson mixing, both in the Standard Model and beyond, are often computed in the 
MS scheme. For convenience, we then also provide here the formulae relating the operators 
renormalized in the RI'-MOM to those renormalized in the MS scheme, at the next-to- 
leading order (i.e. 1-loop). This relation does not depend on the chosen regularization and 
it may be conveniently computed using, for instance, continuum dimensional regularization. 

We restrict our attention to the AF = 2 Parity Conserving operators, which are relevant 
for neutral meson mixing;. These are the 5 operators Qf~ + i of Eq.flU]) with S = +1. 3 For 
these operators, the conversion from the RI'-MOM to the MS scheme can be written in the 
form 

(Q s -- +1 )T=(^^)jQ s -- +1 )™ <«) 



where Ar is a 5 x 5 matrix. This matrix is independent of the choice of the regularization, 
i.e. it is the same for instance for continuum dimensional regularization and for the lattice 
regularization. The chiral symmetry of QCD also implies that the same matrix Ar is also 
valid for the Parity Violating sector (though the operators Q2 and Q2 vanish in the AF = 2 

2 In Ref. indeed, the lattice regularization chosen for the valence quarks s, s', d, d! is maximally twisted 
Wilson fermions with r s = r s < = = — r<j'. In this case, the parity conserving operator Q\ = Ow+aa, 
besides being free from wrong chirality mixings, admits the same renormalization constant (here called 
Zva+av ) as the operator Q\ = Ova+av regularized with untwisted Wilson quarks (i.e. with r s = r s < = 
rd = r d > + 1). 

3 Note that in the AF = 2 case the operators with S = — 1 vanish identically, since 0^ Y = Oxy ■ Moreover, 
there are only three Parity Violating operators, since =+1 and =+1 also vanish. 
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case) . 

When dealing with four-fermion operators, the (modified) minimal subtraction prescrip- 
tion in dimensional regularization is not sufficient however to univocally specify the renor- 
malization scheme. Different MS schemes can be defined, which differ for the definition of 
the so called evanescent operators. The scheme usually adopted in the analysis of K — K 
mixing is the MS scheme defined for instance in Ref. for which the 1-loop conversion 
matrix Ar of Eq. fHHj) reads: 



Ar 



^_M + 8 ln2 










■| - §ln2 -4-41n2 



l-ln2 f-fhi2 

















10 , 10 







1 i 7 



\ 



3 +f ln2 -^ + ^ln2 



56 



3 In 2 -f + fhi2 



(49) 



For B — B mixing, instead, the MS scheme of Ref. [L2| is more commonly adopted. The 
corresponding matrix Ar differ from the one given in Eq. (1491) only in the Qf^ +1 sector, 
which in this case reads 



/ 



Q5 



f + fln2 -^ + ^ln2 



yf + f In 2 -ff + f In 2 



(50) 



In order to correct, to 0(a 2 ), non-perturbative estimates for the renormalization constants 
of 4-fermion operators one should take into account the 0(a 2 ) corrections of Eqs. ( 140]) - (14T|) . 
as well as the 0(a 2 ) terms of the fermion propagator [1]. The exact terms that need to be 
subtracted from the non-perturbative Z q , computed in the RI'-MOM scheme are provided 
in Ref. {3] for general action parameters and in Ref. Q for tree- level Symanzik improved 
gluons, csw = 0, Landau gauge. 
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IV. CONCLUSION 



The calculations presented regard all 4-fermion operators of the form: sXds'Yd' 
where X, Y are generic Dirac matrices. Our results have explicit dependence on: 
p, a, g, csw, A, N c , and implicit depencence on the Symanzik parameters, q. Thus, the 
numerical results are presented for a selection of currently used values of q. 

In a recent paper [4j the present results on the perturbative renormalization of AF = 2 
operators have been combined with numerical simulation data in order to determine non- 
perturbative renormalization coefficients with better precision. This allowed us to extract 
physical values for Bx with reduced lattice artifacts. 



Acknowledgements: Work supported in part by the Research Promotion Foundation 
of Cyprus (Proposal Nr: TEXN/0308/17, ENISX/0506/17). 

Appendix A: Analytic expressions 

In general, the final 1-loop expressions for A^ y , A^ y and A^ y , up to 0(a 2 ), are obtained 
as a function of: the coupling constant g, clover parameter csw-i number of colors N c , lattice 
spacing a, external momentum p, and gauge parameter A. The specific values A = 1 (0) 
correspond to the Feynman (Landau) gauge. Here we present the results for A^ y , A^ Y and 
A^ y for the special choices: r s = = r s > = 7> = 1, and tree-level Symanzik improved 
gluon action. 
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Diagram dl 



\XY(\ 01020304 9 



16tt 2 



a 2 ^ a 



^ x |(A 0(a o ) )^ + a 1 (A e , (ol) )^ + a 2 (A 0(o2) )^}, (Al) 



(A ( a o))d/ = X tl i2 Y l3 u - 0.05294144(3) + 0.737558970(1) c sw + 0.238486988(3) cg w 



-2.100573331(5) A+ -(-1 + A) ln(aV) 



r 2 

-sw 



^(X^)^ J2 (>V) i3 u - 0.507914049(6) + 0.55316919(1) c sw - 0.194516637(3) c 

^(AVVk is (*Wk i 4 [o.018598520(2) - 0.1843897425(8) c sw - 0.0596217473(8) 

1 



// .1' 



— ln(oV 



+ £(x 7 'V) ili2 (y 7 V) i 



0.397715726853 + 0.147715726853 A 



(A2) 



(A 0(al) )^ y = £ (x il , 2 (F 7 ^) isi4 + (A 7 ^) ni2 y i3 U ) x (0.09460083(1)- 0.065711182(4) c sw 



-0.059929106(1) c§w + 0.438508366(3) A 



+ -(-1+ csw-A)ln(aV) 



+ ^ ((AWk l2 (Y^ 4J + (A 7 ") 41 42 (^'Vk u) x (0.1692905881(6) + 0.010283104(5) c sw 



0.0680031615(8) cfw + 0.073857863427 A 



+ — (l + 3csw)ln(oV) 
16 



+ E ((^7 M 7Y), ll2 (W) l3lJ + (^"/Jm^W/)^,) x [ip„(- 0.0279443091(3) csw 



-0.0319830668(5) c| w 

-t^csw ln(a 2 p 2 ) 
16 



(A3) 



17 



(A 0(o 2 ) )J/ = X h l2 Y n u p A (1.32362251(5) - 0.43684285(3) c sw - 0.0208665277(5) c 



sw 



+0.64073441(3) A + ^(-17 + 9 c sw ~ 9 A) In(aV)) 
- ^ 2 Pg (0.06213648(8) - 0.07400055(8) a) 
p 2 (^0.059895142(8) - 0.241755150(3) c sw + 0.114731816(7) cg w 

-0.036928931713 A+ -U-7 + 11 c sw - 4 c| w ) ln(a.'V)) 
7^(1.01694823(2) - 0.44474062(1) c sw - 0.033265121(3) Cg W 



I- 1 



(o.00592406(2) - 0.00295805(2) . 



^(*7 M k * 2 (*Vk u [p»Pv ( - 0.19915360(1) + 0.212823513(3) c sw + 0.033028338(2) c 



2 

SW 



+0.1600141922(8) A + ^(-4 + 5 c sw + 2 c^ w - 3 A) ln(a V 



^(-XWk l2 (l^y'Vk 14 [p 2 ( - 0.08962805(1) + 0.0769373498(3) c sw + 0.0067184623(3) Cg W 



-_(7-5c S w)ln(aV) 



p 2 ( + 0.16608907(6) + 0.07446360(2) c sw - 0.0087763322(3) c 



sw 



29 
^180 

o 



ln(aV)) 



- 0.048180849735 



^ (X 7 'V) il t2 (F 7 Vk i4 2VV ( - 0.21865904(4) + 0.054629909(8) c sw + 0.00276900638(7) c 

-0.082411837(6) A + ^7j( 164 _ 60 c sw + 45 A) ln(aV)) 



{ 



0.110138789528 - 0.024619287809 A 



) 



+P^ TW (0.140961390102 + 0.045240352404 A 



P ^r R ( _ 0.477634781(8) - 0.083831642(8) a) 
(Wk * 2 ^ [p* (0.00385492 



£ (jfyYYk i2 (y 7 "yy ) ls u ^ Pv ( - 0.0209503296(6) c 2 w - 1 c 2 w In(a V)) ■ 



Diagram d2 



AfW^Sr = ^5a ia Ja 3 a 4 (n c - x {(A 0(a o))f/ + a\k ^))l Y + a 2 (A 0(q2) )l Y I , (A5) 



(Ao(aO))f 2 F 



X 4l i2 F i3 14 1.2904478(4) + 0.737558970(1) c sw + 0.238486988(3) c| w 



+2.3960046(4) A + -(-1 - A) ln(aV) 



A n i2 (7^7^., l4 \ - 0.507914047(8) + 0.55316917(2) c sw - 0.194516638(9) c| w 
^ X 4l i2 (YV*W)i 3 i 4 [ ^ 0.129117207(2) - 0.1843897425(8) c sw - 0.0596217473(8) c§ w 



//.// 



-- ln(a 2 p 2 



+ ^x JlJ2 (7^y 7 V) ! : 





-->)} 


I p 2 v 


4 /J 



(A6) 



(A 0(Ql) )^ = £ (x lll2 (r 7 ^) 43? ; 4 +x llt2 (7^r) 4 



(0.37785613(9) - 0.56675680(2) 



csw 



-0.160026205(2) cj w - 0.48393977(9) A 



--(-3 + 5 c sw + 2 A) ln(aV; 



+ Y, fci !2 (7Y5 / 7 i/ k ! 4+^M 2 (fW)M 4 ) x 0.073251555(1) -0.001704761(4) c sw 



-0.032093938(2) c 



sw 



+ 77(3- csw)ln(aV) 
16 



E (An J2 (7 M 7 l/ 7^7"7 P ),3»4+A ni2 (7 l/ 7 P y7 l/ 7V) J3 ,4) x [i Pfl (o.0459135542(5) c sw 



+0.0319830668(5) c| w 
1 



c S w ln(a 2 p 2 



16 



(A7) 



19 



(^o(ai))d 2 Y = X h l2 Y l3 u [p 2 (o.7374671(6) - 0.24301094(4) c sw - 0.0096054476(8) 



; sw 



-0.4696085(6) A + —(-23 + 5 



(^- 77+15A ))] 



c sw + 15 A) ln(aV)) 



- 51 l 2 (7 M *Vk i4 P 2 (0.04610701(4) - 0.19136171(3) c sw + 0.02347831(4) 



'SW 



-0.06249999(1) A + —(1 + c sw )ln(aV)) 
16 / 

+pl (o.17251518(3) - 0.19211806(3) c sw + 0.01744902(2) cfw) 



i 2 i 



(lWku+* ; 
^X ll . l2 (7^7 1 ') 



•P^/101 



- P 



/101N 

V288/ 



i 3 i 4 



p M p„ ^0.05513763(3)- 
-0.10887203(3) ' 



0.005630284(8) c sw - 0.072690409(7) 
A + ^(-2 + c sw + c| w + A) ln(aV) 



-sw 



? ( - 0.05064893(1) + 0.0394316274(6) c sw + 0.00332360968(9) 

+ ^(13-15c sw )ln(aV)) 
p 2 ( + 0.05383442(9) + 0.124311493(7) c sw - 0.0051958638(1) c 

4 h <" v >) 

V 240/. 



c sw 



x n h il >1 l PY l P l v )i3 u V^Vv ( - 0.03270359(5) - 0.039026988(4) c sw + 0.0015068706(2) - 

-0.06926696(4) A + y^(28 + 60 c sw + 45 A) ln(«V)) 



'SW 



£ (0.00385492^5(2) c 2 w ) 

(0.015978597(1) c| w - ^ c 2 w ln(a 2 p 2 )) 



£ X n^(7 



- E ^ 



Diagram d3 



16tt 2 



a^a 3 a2 



a2 ^(13 a4 



x (A 0(a o ) )f 3 y + a 1 (A 0(a i ) )f 3 r + a 2 (A 0(Q 2 ) )^ 



(A9) 



(Ao(qO) 



1.2904478(4) + 0.737558970(1) c sw + 0.238486988(3) c| w 



+2.3960046(4) A + -(-1 - A) ln(aV) 



.2 

-sw 



^(X 7 ^) 4l i2 (7 M F) l3 u - 0.507914047(8) + 0.55316917(2) c sw - 0.194516638(9) c; 
X)(*W)ii i 2 (l>*YY) l3 44 [ - 0.129117207(2) - 0.1843897425(8) c sw - 0.0596217473(8) c 



,2 

SW 



/ / . v 



+ iln(aV) 





-->)] 


l p 2 V 


4 /J 



(A10) 



1 c S w 



(A ( fl i))f =E ((^7 M )iii 2 ^3u +X ili2 (^Y) i3ii ) x [ip„ (0.37785613(9) -0.56675680(2), 

— 0.160026205(2) C g W - 0.48393977(9) A 
+ i(-3 + 5c sw + 2A) ln(aV)) 

+ H (( x ^hi^2(^YYh 3 ll + (X^^) lll2 {^Y) l3l ^ x [ip M (- 0.073251555(1) -0.001704761(4) c sw 



//./-' 



-0.032093938(2) c 



sw 



+ttt(3- csw)ln(aV) 
16 



+ E Y/),: » 2 (f 7^),, I4 + (^7 Y).: !2 (7"/7^) t 3 u) x (0.0459135542(5) c sw 



-0.0319830668(5) c£w 

-777 c S w ln(aV 
16 



(AH) 



21 



(A (a 2 ))f 3 y = X n h Y i 3 u P 2 (o. 7374671(6) - 0.24301094(4) c sw - 0.0096054476(8) 

-0.4696085(6) A + ^^(-23 + 5 c sw + 15 A) In(aV) 



c sw 



120 

"^(^(-77+ 15 A) 



^ * 2 (7 M ^k ,4 [p 2 (0.04610701(4) - 0.19136171(3) c sw + 0.02347831(4) c 2 sw 

— 0.06249999(1) A + — (1 + c sw ) ln(aV) 
16 



c sw 



+p^(o. 17251518(3) - 0.19211806(3) c sw + 0.01744902(2) 

^((x 7 Y), 1 , ! i',3, 4 +i, 1 , ! (TYi'k i4 ) x 

- Y^iX^K *2 (l v Y)i 3 u [p/iP^ (0.05513763(3) - 0.005630284(8) c sw - 0.072690409(7) 



-sw 



-0.10887203(3) A + ^(—2 + c sw + c| w + A) In(aV) 
8 



■^{XjU'fh ^ (j^YY) h u \p 2 ( - 0.05064893(1) + 0.0394316274(6) c sw + 0.00332360968(9) c§ w 

+^(13-15c sw )ln(aV)) 



( + 0.05383442(9) + 0.124311493(7) c sw - 0.0051958638(1) c| w 



1 



5 In(oV) 
240 



P 



{XYl^n ^2 (l"Y Y ) la u P^Pv ( - 0.03270359(5) - 0.039026988(4) c sw + 0.0015068706(2) c 



2 

SW 



-0.06926696(4) A + ^^(28 + 60 c sw + 45 A) ln(aV 



(pIpv+PhpI) ( 1 



f (41 -40 A) 

V960 V ' 



P»PV ^2)2 



(%5< 7 + 25A > 



^(i<-*-^)) 

£ (XfVVk i a (7 M 7V30ia .4 [^(0.0038549^95(2) c^ w 
^ (X 7 ^ 7 " 7 CT ) ll l2 ( 7 " 7 V V) i3 » 4 [p^„ (0.015978597(1) c§ w - ^ c§ w ln(a V) 



(A12) 



Appendix B: Notation in ASCII file: 4-fermi.m 



The full body of our results can be accessed online through the file 4-fermi.m, which is a 
Mathematica input file. It includes the expressions for the three Feynman diagrams: 

• A^ Y : dl [action, csw,lambda,Nc,g,aL] 

• A^ y : d2 [action, csw,lambda,Nc,g,aL] 

• A^ y : d3 [action, csw,lambda,Nc,g,aL] 

from which one can construct the matrix elements of any 4-fermion operator of the above 
form. Each expression depends on the variables: 

• action: Selection of improved gauge action as follows, 1 — > Plaquette, 2 — > Tree Level 
Symanzik, 3 ->■ TILW (f3 c = 8.60), 4 -> TILW (f3 c = 8.45), 5 ->■ TILW (f3 c = 8.30), 
6 ->• TILW (pc = 8.20), 7 ->■ TILW (/3c = 8.10), 8 ^ TILW (f3 c = 8.00), 9 ->■ 
Iwasaki, 10 -> DBW2 

• csw: clover parameter 

• lambda: gauge parameter (Landau/Feynman/Generic correspond to 0/1/lambda) 

• Nc: number of colors 

• g: coupling constant 

• aL: lattice spacing 

In particular, the quantities of interest in that file are the renornalization matrices for the 
10 Parity Conserving operators, and 10 Parity Violating operators, which read: 

• Z s=+1 : PCplus[action,csw,lambda,Nc,g,aL,p2,p4] [Projector, LGreen] 

• Z s= ~ x : PCminusfaction, csw, lambda, Nc,g,aL,p2,p4] [Projector, LGreen] 

• Z s=+1 : PVplus[action,csw,lambda,Nc,g,aL,p2,p4] [Projector, LGreen] 

• Z s= ~ x : PVminus [action, csw, lambda,Nc,g,aL,p2,p4] [Projector, LGreen] . 

The additinal variables are 

• P2: EtiP? 
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p 4: EtiPi 

Projector: the index I of Section HTT1 (1 to 5) 
LGreen: the index m of Section HTT1 (1 to 5) 
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Parity Conserving, Flavour Exchange Symmetry Plus 
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2 


-8/3 


-2.79899092(5) 


0.39295395(5) 


(1, 


2) 


5.41774985(8) 


-5.9004712(3) 


2.0748441(1) 

















(1, 


3) 


-0.67721873(1) 


0.73755889(5) 


-0.25935552(2) 

















(1, 


4) 


-0.67721873(1) 


0.73755892(2) 


-0.259355516(6) 

















(1, 


5) 


-2.03165620(4) 


2.21267677(7) 


-0.77806655(2) 

















(2, 


1) 


7.4494060(1) 


-8.1131478(4) 


2.8529107(1) 

















(2, 


2) 


4.331175(2) 


-3.68779471(5) 


0.30109257(1) 


10.816593(2) 


1 


-8/3 


-0.87361423(2) 


0.57016434(2) 


(2, 


3) 


-0.8545378(4) 


-0.368779461(8) 


0.129677758(2) 


-0.1931470(4) 








0.89270365(2) 


0.08962151(2) 


(2, 


4) 


0.338609366(5) 


-0.36877945(2) 


0.129677758(6) 

















(2, 


5) 


-1.01582810(1) 


1.10633834(5) 


-0.38903328(2) 

















(3, 


1) 


1.35443746(2) 


-1.47511779(6) 


0.51871103(2) 

















(3, 


2) 


-9.615778(1) 


-10.32582548(4) 


-2.34313283(3) 


-3.772588(1) 


6 





1.70275904(9) 


0.92098603(9) 


(3, 


3) 


13.627614(2) 


9.58826675(6) 


4.59385837(5) 


15.316593(2) 


-8 


-8/3 


1.92846910(2) 


-0.27358566(2) 


(3, 


4) 


-8.8038435(2) 


9.5882656(6) 


-3.3716217(2) 

















(3, 


5) 


-6.09496858(8) 


6.6380301(3) 


-2.3341997(1) 

















(4, 


1) 


-2.70887493(5) 


2.9502357(1) 


-1.03742206(3) 

















(4, 


2) 


9.4810622(1) 


-10.3258245(4) 


3.6309772(1) 

















(4, 


3) 


-10.8354997(2) 


11.8009423(6) 


-4.1496883(2) 

















(4, 


4) 


10.269733(2) 


7.37558970(4) 


2.38486989(5) 


14.816594(2) 


-5 


-8/3 


2.01567488(3) 


0.02849782(3) 


(4, 


5) 


9.732709(2) 


7.37558970(3) 


2.38486988(5) 


3.477157(2) 


-5 





-1.07855468(7) 


-0.76054387(7) 


(5, 


1) 


-2.70887493(5) 


2.9502357(1) 


-1.03742206(2) 

















(5, 


2) 


-1.35443746(2) 


1.47511779(6) 


-0.51871103(2) 

















(5, 


3) 


























(5, 


4) 


1.1783609(6) 


-0.49170598(2) 


-0.15899133(2) 


0.1590523(6) 


1/3 





-0.98220341(2) 


-0.06601462(2) 


(5, 


5) 


2.297078(2) 


-8.35900166(4) 


-2.70285255(5) 


10.134698(2) 


17/3 


-8/3 


-3.06215787(3) 


0.83396857(3) 



TABLE I: The coefficients - <i°' 6) and ~ 4™. 



Parity Conserving, Flavour Exchange Symmetry Plus 



(I, m) 


,+ (2,3) 
a l,m 


r+(2,4) 

a l,m 


,+("2,51 
a l,m 


,+(2,61 
a l,m 


.+("2,71 
a l,m 


,+(2,8) 


,+(2,9) 
a l,m 


,+ (2 10) 
u l,m 


(1, 1) 


2.642227(4) 


1.7473718(2) 


0.08346609(1) 


-3.289500(4) 


-19/18 


-1/2 





25/24 


(1, 2) 


-1.0988822(5) 


2.5685206(3) 


-0.9493554(4) 


0.9569921(2) 





-1 


1 


-1/3 


(!> 3) 


0.13736030(9) 


-0.32106507(4) 


0.11866942(8) 


-0.11962402(4) 





1/8 


-1/8 


1/24 


(1, 4) 


0.35245840(2) 


-0.399645902(8) 


0.06172133(2) 


0.0040994899(4) 


-1/4 


3/8 


-1/4 


-1/24 


(!> 5) 


1.05737519(6) 


-1.19893771(2) 


0.18516399(5) 


0.012298470(1) 


-3/4 


9/8 


-3/4 


-1/8 


(2, 1) 


-1.5109631(5) 


3.5317158(3) 


-1.3053636(4) 


1.3158642(3) 





-11/8 


11/8 


-11/24 


(2, 2) 


0.789419(4) 


1.84106755(5) 


-0.37712545(3) 


-2.739903(4) 


1/18 


-11/8 


3/8 


41/48 


(2, 3) 


-1.6447718(6) 


0.199822951(3) 


-0.065051477(5) 


-0.0214962(6) 


11/24 


-3/16 


1/16 


-1/24 


(2, 4) 


-0.06868015(3) 


0.16053253(1) 


-0.01966766(2) 


0.05981201(1) 





-1/16 





-1/48 


(2, 5) 


0.20604045(9) 


-0.48159760(4) 


0.05900298(7) 


-0.17943602(4) 





3/16 





1/16 


(3, 1) 


-0.2747206(1) 


0.64213013(6) 


-0.2373388(1) 


0.23924803(5) 





-1/4 


1/4 


-1/12 


(3, 2) 


-7.416224(3) 


4.6520712(2) 


-0.18110545(2) 


0.870210(3) 


5/3 


-9/4 


1/4 


-13/24 


(3, 3) 


2.045124(4) 


-3.9381018(3) 


-0.49577614(3) 


-4.174196(4) 


11/36 


7/8 


3/8 


17/12 


(3, 4) 


1.785683(1) 


-4.1738459(5) 


0.5113591(7) 


-1.5551122(4) 





13/8 





13/24 


(3, 5) 


1.2362427(5) 


-2.8895856(3) 


0.3540179(4) 


-1.0766161(2) 





9/8 





3/8 


(4, 1) 


1.40983359(8) 


-1.59858361(4) 


0.24688532(7) 


0.016397959(2) 


-1 


3/2 


-1 


-1/6 


(4, 2) 


-1.9230442(7) 


4.4949109(3) 


-0.5506945(6) 


1.6747362(3) 





-7/4 





-7/12 


(4 5 3) 


2.197764(1) 


-5.1370411(5) 


0.6293651(7) 


-1.9139843(4) 





2 





2/3 


(4, 4) 


-0.286209(4) 


-2.8633160(3) 


-0.023092586(8) 


-4.273362(4) 


1/18 


5/4 





59/48 


( 4 > 5 ) 


4.602710(3) 


-3.0948716(3) 


-0.05164224(1) 


-0.928529(3) 


-7/9 


5/4 





9/16 


(5, 1) 


1.40983359(8) 


-1.59858361(3) 


0.24688532(6) 


0.016397959(2) 


-1 


3/2 


-1 


-1/6 


(5, 2) 


0.2747206(1) 


-0.64213013(6) 


0.0786706(1) 


-0.23924803(5) 





1/4 





1/12 


(5, 3) 


























(5, 4) 


1.255191(1) 


0.2526361(1) 


0.009152736(3) 


0.009222(1) 


-17/54 


-1/12 





1/16 


(5, 5) 


0.828945(4) 


4.0632560(3) 


0.12704697(1) 


-2.661259(4) 


-23/27 


-17/12 





35/48 



TABLE II: The coefficients d^' 3) - d} 



Parity Conserving, Flavour Exchange Symmetry Minus 



(I, m) 


,-(0,1) 
a l,m 


i— (0,2) 
a l,m 


,-(0,3) 
a l,m 


,-(0,4) 
a l,m 


,-(0,5) 
A ^ ' > 

a l,m 


,-(0,6) 
a l,m. 


,-(2 1) 
a l,m 


,-(2,2) 
a l,m 


(1, 1) 


-2.830716(3) 


5.90047176(3) 


1.90789592(4) 


9.748573(3) 


-4 


-8/3 


2.12575961(5) 


0.46409210(5) 


(1, 2) 


5.41774985(8) 


-5.9004712(3) 


2.0748441(1) 

















(!> 3) 


1.35443746(2) 


-1.47511779(6) 


0.51871103(2) 

















(1, 4) 


1.35443746(2) 


-1.47511784(3) 


0.518711032(7) 

















(!> 5) 


4.06331239(5) 


-4.42535353(9) 


1.55613310(2) 

















(2, 1) 


3.38609366(9) 


-3.6877945(3) 


1.2967776(1) 

















(2, 2) 


0.267862(2) 


0.73755883(5) 


-1.25504053(1) 


10.816593(2) 


1 


-8/3 


-0.87361423(2) 


0.57016434(2) 


(2, 3) 


1.5317566(4) 


-0.368779461(8) 


0.129677758(2) 


0.1931470(4) 








-0.89270365(2) 


-0.08962151(2) 


(2, 4) 


0.338609366(5) 


-0.36877945(2) 


0.129677758(6) 

















(2, 5) 


-1.01582810(1) 


1.10633834(5) 


-0.38903328(2) 

















(3, 1) 


1.35443746(2) 


-1.47511779(6) 


0.51871103(2) 

















(3, 2) 


12.324653(1) 


7.37558979(4) 


3.38055490(3) 


3.772588(1) 


-6 





-1.70275904(9) 


-0.92098603(9) 


(3, 3) 


9.564301(2) 


14.01362029(6) 


3.03772527(5) 


15.316593(2) 


-8 


-8/3 


1.92846910(2) 


-0.27358566(2) 


(3, 4) 


-12.8671559(2) 


14.0136190(6) 


-4.9277548(2) 

















(3, 5) 


6.09496858(8) 


-6.6380301(3) 


2.3341997(1) 

















(4, 1) 


5.41774986(5) 


-5.9004714(1) 


2.07484413(3) 

















(4, 2) 


-6.7721873(1) 


7.3755889(4) 


-2.5935552(1) 

















(4 5 3) 


-10.8354997(2) 


11.8009423(6) 


-4.1496883(2) 

















(4, 4) 


12.922182(2) 


16.22629734(4) 


5.24671376(5) 


15.816593(2) 


-11 


-8/3 


1.84126332(3) 


-0.57566914(3) 


( 4 > 5 ) 


-9.555272(2) 


-1.47511794(3) 


-0.47697397(5) 


-3.068020(2) 


1 





2.15255184(7) 


0.47726124(7) 


(5, 1) 


5.41774986(6) 


-5.9004714(1) 


2.07484413(3) 

















(5, 2) 


-1.35443746(2) 


1.47511779(6) 


-0.51871103(2) 

















(5, 3) 


























(5, 4) 


-1.1192154(6) 


2.45852990(2) 


0.79495663(2) 


-0.0226732(6) 


-5/3 





1.34020247(2) 


-0.02841292(2) 


(5, 5) 


-3.777376(2) 


0.49170598(4) 


0.15899132(5) 


8.771247(2) 


-1/3 


-8/3 


0.16287196(3) 


0.68000502(3) 



TABLE III: The coefficients d^ 0,1) - d, ^°' 6) and - 



Parity Conserving, Flavour Exchange Symmetry Minus 



(I, m) 


i— (2,31 


,-(2,4) 
a l,m 


,—(2,51 
a l,m 


,— (2,61 
a l,m 


,-(2,7) 
a l,m 


,-(2,8) 
a l,m 


,-(2 9) 
a l,m 


,-(2 10) 
a l,m 


(1, 1) 


3.611581(4) 


-3.4947437(3) 


-0.16693217(1) 


-2.043017(4) 


-5/9 


1 





5/12 


(1, 2) 


-1.0988822(5) 


2.5685206(3) 


-0.9493554(4) 


0.9569921(2) 





-1 


1 


-1/3 


(!> 3) 


-0.2747206(1) 


0.64213013(6) 


-0.2373388(1) 


0.23924803(5) 





-1/4 


1/4 


-1/12 


(1, 4) 


-0.70491680(2) 


0.79929180(1) 


-0.12344266(2) 


-0.0081989797(6) 


1/2 


-3/4 


1/2 


1/12 


(!> 5) 


-2.11475039(7) 


2.39787541(3) 


-0.37032798(6) 


-0.024596939(2) 


3/2 


-9/4 


3/2 


1/4 


(2, 1) 


-0.6868013(5) 


1.6053254(3) 


-0.5933471(4) 


0.5981201(2) 





-5/8 


5/8 


-5/24 


(2, 2) 


2.904169(4) 


-0.55680786(5) 


0.40349227(3) 


-2.715306(4) 


-13/9 


7/8 


-3/8 


29/48 


(2, 3) 


1.2923134(6) 


0.199822951(3) 


-0.065051477(5) 


0.0173967(6) 


-5/24 


-3/16 


1/16 


1/12 


(2, 4) 


-0.06868015(3) 


0.16053253(1) 


-0.01966766(2) 


0.05981201(1) 





-1/16 





-1/48 


(2, 5) 


0.20604045(9) 


-0.48159760(4) 


0.05900298(7) 


-0.17943602(4) 





3/16 





1/16 


(3, 1) 


-0.2747206(1) 


0.64213013(6) 


-0.2373388(1) 


0.23924803(5) 





-1/4 


1/4 


-1/12 


(3, 2) 


6.006391(3) 


-3.0534875(2) 


-0.33930637(2) 


-0.886608(3) 


-2/3 


3/4 


1/4 


17/24 


(3, 3) 


4.159875(4) 


-6.3359772(3) 


0.28484158(3) 


-4.149599(4) 


-43/36 


25/8 


-3/8 


7/6 


(3, 4) 


2.609845(1) 


-6.1002363(5) 


0.7473711(7) 


-2.2728563(4) 





19/8 





19/24 


(3, 5) 


-1.2362427(5) 


2.8895856(3) 


-0.3540179(4) 


1.0766161(2) 





-9/8 





-3/8 


(4, 1) 


-2.81966719(8) 


3.19716722(4) 


-0.49377065(7) 


-0.032795919(2) 


2 


-3 


2 


1/3 


(4, 2) 


1.3736030(7) 


-3.2106507(3) 


0.3933532(6) 


-1.1962402(3) 





5/4 





5/12 


(4 5 3) 


2.197764(1) 


-5.1370411(5) 


0.6293651(7) 


-1.9139843(4) 





2 





2/3 


(4, 4) 


6.491207(4) 


-7.4107631(3) 


-0.187841964(8) 


-4.050432(4) 


-17/18 


11/4 





65/48 


( 4 > 5 ) 


-2.042489(3) 


0.0632400(3) 


-0.05819067(1) 


1.051220(3) 


5/9 


-1/4 





-9/16 


(5, 1) 


-2.8196672(1) 


3.19716722(4) 


-0.49377065(8) 


-0.032795919(2) 


2 


-3 


2 


1/3 


(5, 2) 


0.2747206(1) 


-0.64213013(6) 


0.0786706(1) 


-0.23924803(5) 





1/4 





1/12 


(5, 3) 


























(5, 4) 


-0.401784(1) 


-1.2631800(1) 


-0.045763706(3) 


0.031675(1) 


13/54 


5/12 





-1/16 


(5, 5) 


3.734320(4) 


-0.9473033(3) 


-0.09480166(1) 


-1.755961(4) 


-23/27 


1/12 





17/48 



TABLE IV: The coefficients d^ 2,3) - dj 



Parity Violating, Flavour Exchange Symmetry Plus 



(I, m) 


r+(0,l) 
l,m 


£+(0,2) 
°l,m 


£+(0,3) 
l,m 


£+(0,4) 
l,m 


^+(0,5) 
°l,m 


£+(0,6) 
°l,m 


£+(2,D 
i,m 


£+(2,2) 
l,m 


(1, 1) 


7.607190(2) 


-2.95023588(3) 


-0.95394796(3) 


12.293750(2) 


2 


-8/3 


-2.79899092(5) 


0.39295395(5) 


(2, 2) 
(2, 3) 


2.299519(2) 
-1.1931472(4) 


-1.475117940(5) 



-0.476973978(5) 



10.816593(2) 
-0.1931470(4) 


1 




-8/3 



-0.87361423(2) 
0.89270365(2) 


0.57016434(2) 
0.08962151(2) 


(3, 2) 
(3, 3) 


-10.970215(1) 
11.595958(2) 


-8.85070764(3) 
11.80094352(4) 


-2.86184387(3) 
3.81579182(4) 


-3.772588(1) 
15.316593(2) 


6 
-8 




-8/3 


1.70275904(9) 
1.92846910(2) 


0.92098603(9) 
-0.27358566(2) 


(4, 4) 
(4, 5) 


10.269733(2) 
9.732709(2) 


7.37558970(4) 
7.37558970(3) 


2.38486989(5) 
2.38486988(5) 


14.816594(2) 
3.477157(2) 


-5 
-5 


-8/3 



2.01567488(3) 
-1.07855468(7) 


0.02849782(3) 
-0.76054387(7) 


(5, 4) 
(5, 5) 


1.1783609(6) 
2.297078(2) 


-0.49170598(2) 
-8.35900166(4) 


-0.15899133(2) 
-2.70285255(5) 


0.1590523(6) 
10.134698(2) 


1/3 
17/3 




-8/3 


-0.98220341(2) 
-3.06215787(3) 


-0.06601462(2) 
0.83396857(3) 



TABLE V: The coefficients stf* - and 6^ - 5^. 



Parity Violating, Flavour Exchange Symmetry Plus 



(I, m) 


£+(2,3) 
°l,m 


£+(2,4) 
°l,m 


£+(2,5) 
°l,m 


£+(2,6) 
°l,m 


r+(2,7) 
°l,m 


A +(2,8) 
°l,m 


£+(2,9) 
°l,rn 


£+(2,10) 
°l,m 


(1, 1) 


2.642227(4) 


1.7473718(2) 


0.08346609(1) 


-3.289500(4) 


-19/18 


-1/2 





25/24 


(2, 2) 
(2, 3) 


1.846794(4) 
-1.4685426(6) 


0.64212984(5) 



0.0131834100(7) 



-2.727604(4) 
-0.0194464(6) 


-25/36 
1/3 


-1/4 







35/48 
-1/16 


(3, 2) 
(3, 3) 


-6.711307(3) 
3.102499(4) 


3.8527794(2) 
-5.1370395(3) 


0.079100460(4) 
-0.105467280(7) 


0.878409(3) 
-4.161897(4) 


7/6 
-4/9 


-3/2 
2 






-5/8 
31/24 


(4, 4) 
(4, 5) 


-0.286209(4) 
4.602710(3) 


-2.8633160(3) 
-3.0948716(3) 


-0.023092586(8) 
-0.05164224(1) 


-4.273362(4) 
-0.928529(3) 


1/18 
-7/9 


5/4 
5/4 






59/48 
9/16 


(5, 4) 
(5, 5) 


1.255191(1) 
0.828945(4) 


0.2526361(1) 
4.0632560(3) 


0.009152736(3) 
0.12704697(1) 


0.009222(1) 
-2.661259(4) 


-17/54 
-23/27 


-1/12 
-17/12 






1/16 

35/48 



TABLE VI: The coefficients 5^ ,3) - 5] 



Parity Violating, Flavour Exchange Symmetry Minus 



(I, m) 


,-(o,i) 


,-(0,2) 


,-(0,3) 


,-(o,4) 


,-(o,5) 

°l,m 


,-(0,6) 
°l,m 


57™ 

l,m 


,-(2,2) 


(1, 1) 


-2.830716(3) 


5.90047176(3) 


1.90789592(4) 


9.748573(3) 


-4 


-8/3 


2.12575961(5) 


0.46409210(5) 


(2, 2) 
(2, 3) 


2.299519(2) 
1.1931472(4) 


-1.475117940(5) 



-0.476973978(5) 



10.816593(2) 
0.1931470(4) 


1 




-8/3 



-0.87361423(2) 
-0.89270365(2) 


0.57016434(2) 
-0.08962151(2) 


(3, 2) 
(3, 3) 


10.970215(1) 
11.595958(2) 


8.85070764(3) 
11.80094352(4) 


2.86184387(3) 
3.81579182(4) 


3.772588(1) 
15.316593(2) 


-6 
-8 




-8/3 


-1.70275904(9) 
1.92846910(2) 


-0.92098603(9) 
-0.27358566(2) 


(4, 4) 
(4, 5) 


12.922182(2) 
-9.555272(2) 


16.22629734(4) 
-1.47511794(3) 


5.24671376(5) 
-0.47697397(5) 


15.816593(2) 
-3.068020(2) 


-11 
1 


-8/3 



1.84126332(3) 
2.15255184(7) 


-0.57566914(3) 
0.47726124(7) 


(5, 4) 
(5, 5) 


-1.1192154(6) 
-3.777376(2) 


2.45852990(2) 
0.49170598(4) 


0.79495663(2) 
0.15899132(5) 


-0.0226732(6) 
8.771247(2) 


-5/3 
-1/3 




-8/3 


1.34020247(2) 
0.16287196(3) 


-0.02841292(2) 
0.68000502(3) 



TABLE VII: The coefficients 6~™ - 5~^ ] and - 5~ 



Parity Violating, Flavour Exchange Symmetry Minus 



(I, m) 


,-(2,3) 
°l,m 


,-(2,4) 
°l,m 


,-(2,5) 
°l,m. 


,-(2,6) 
°l,m. 


,-(2,7) 

°l,m. 


,-(2,8) 
°l,m. 


,-(2,9) 
°l,rn 


,-(2,10) 
°l,m 


(1, 1) 


3.611581(4) 


-3.4947437(3) 


-0.16693217(1) 


-2.043017(4) 


-5/9 


1 





5/12 


(2, 2) 


1.846794(4) 


0.64212984(5) 


0.0131834100(7) 


-2.727604(4) 


-25/36 


-1/4 





35/48 


(2, 3) 


1.4685426(6) 








0.0194464(6) 


-1/3 








1/16 


(3, 2) 


6.711307(3) 


-3.8527794(2) 


-0.079100460(4) 


-0.878409(3) 


-7/6 


3/2 





5/8 


(3, 3) 


3.102499(4) 


-5.1370395(3) 


-0.105467280(7) 


-4.161897(4) 


-4/9 


2 





31/24 


(4, 4) 


6.491207(4) 


-7.4107631(3) 


-0.187841964(8) 


-4.050432(4) 


-17/18 


11/4 





65/48 


(4, 5) 


-2.042489(3) 


0.0632400(3) 


-0.05819067(1) 


1.051220(3) 


5/9 


-1/4 





-9/16 


(5, 4) 


-0.401784(1) 


-1.2631800(1) 


-0.045763706(3) 


0.031675(1) 


13/54 


5/12 





-1/16 


(5, 5) 


3.734320(4) 


-0.9473033(3) 


-0.09480166(1) 


-1.755961(4) 


-23/27 


1/12 





17/48 



TABLE VIII: The coefficients - S t j 



